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IFD:  WIDE-ANGLE  CAPABILITY 

1.  INTRODUCTION 


A  useful  approach  for  calculating  underwater-sound  propagation  in  a 
range-dependent  environment  was  introduced  by  Tappert,!-^  who  used  a  para¬ 
bolic-equation  (PE)  method  to  solve  the  acoustic  wave  equation.  Tappert's  PE 
approximation  transforms  the  acoustic  wave  equation,  which  is  an  elliptic 
partial-differential  equation,  into  a  parabolic  partial -differential  equation. 
The  resulting  PE  governs  the  outgoing  field  and  can  be  solved  by  the  split-step 
Fourier  algorithm  developed  by  Tappert.  The  split-step  algorithm  solves  the 
parabolic  wave  equation  by  imposing  an  artificial  zero  bottom-boundary  condi¬ 
tion  and  pressure-release  surface  condition;  therefore,  it  solves  a  pure  ini¬ 
tial-value  problem.  In  cases  where  bottom  interaction  is  strong,  a  more 
general-purpose  solution  is  useful  for  handling  a  bottom-boundary  condition. 

For  this  reason,  Lee,  Botseas,  and  Papadakis8>°  developed  an  implicit  finite- 
difference  (IFD)  scheme  as  well  as  an  ordinary-differential  (ODE)  method  to 
solve  Tappert's  parabolic  wave  equation.  Their  schemes  are  designed  to  treat 
the  bottom-boundary  condition  exactly.  Carnahan's?  approach  in  handling  the 
horizontal-interface  boundary  condition  for  the  heat-transfer  problem  was  used 
by  McDaniel  and  Lee8  to  incorporate  interface  conditions  into  the  parabolic 
wave  equation. 

These  interface  conditions  were  programmed  into  the  IFD  computer  model 
by  Lee  and  Botseas.9  Since  then,  GilbertlO  and  Greenell  derived  an  improved 
equation  that  extends  the  treatment  of  propagation  angles  in  Tappert's  equa¬ 
tion  from  less  than  approximately  15  deg  up  to  approximately  40  deg.  For 
wide-angle  (>15  deg)  propagation,  Tappert's  original  equation  can  introduce 
large  phase  errors  into  the  solution;  furthermore,  it  is  with  wide-angle  prop¬ 
agation  that  strong  bottom  interaction  generally  takes  place.  Consequently, 
accurate  treatment  of  bottom-boundary  conditions  and  wide-angle  propagation 
are  closely  related  problems.  Recent  progress  in  the  theory,  development,  and 
implementation  of  a  wide-angle  PE  is  reported  in  references  12  and  15.  The 
incorporation  of  the  wide-angle  capability  into  the  IFD  model  is  the  topic  of 
this  report. 

The  mathematical  formulation  and  modifications  that  incorporate  the  wide- 
angle  capability  into  the  IFD  model  are  given  in  the  next  section.  The  same 
approach  used  by  McDaniel  and  Lee8  is  used  to  incorporate  the  interface  con¬ 
ditions  into  the  improved  wide-angle  parabolic  equation.  Section  3  describes 
input  and  output  formats  in  detail,  while  section  4  is  devoted  to  a  test  prob¬ 
lem. 


External  users  can  obtain  a  copy  of  the  improved  model  by  forwarding  a 
formal  request  and  a  blank  magnetic  tape  to  the  Naval  Underwater  Systems  Cen¬ 
ter  (NUSC) .  All  programs  are  written  in  FORTRAN  for  the  VAX-11/780  computer. 

It  is  requested  that  the  authors  be  notified  if  any  difficulties  with  the 
model  are  experienced.  User  contributions  that  will  enhance  the  model  are 
invited . 
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I.  MATHEMATICAL  FORMULATION 


The  derivation  of  the  improved  wide-angle  parabolic  equation  is  presented 
in  four  parts.  First,  McDaniel  and  Lee's8  discussion  on  horizontal-interface 
conditions  that  apply  to  the  two-way  (outgoing  and  incoming)  farfield  equation 
is  presented.  Second,  McDaniel  and  Lee's8  procedure  is  used  to  incorporate 
these  interface  conditions  into  the  two-way  farfield  equation  for  the  acoustic 
field.  Third,  the  farfield  equation  with  interface  conditions  is  formally 
solved,  resulting  in  an  improved  one-way  outgoing  wave  equation.  Finally,  the 
coefficients  that  give  the  improved  equation  wide-angle  capability  are  pre¬ 
sented. 


2.1  HORIZONTAL- INTERFACE  CONDITIONS 

Consider  two  different  media,  as  shown  in  figure  1,  in  which  a  numerical 
subscript  indicates  the  medium,  and  is  the  i,  tort  ice  depth.  Let  o  be  a 
constant  density,  and  let  p  be  the  acoustic  pressure,  lor  a  cy 1 indrical ly 
symmetric  geometry  and  a  harmonic  source,  p  sat i sties  the  reduced  wave  equa¬ 
tion 

72p  *  k'jn'l.r.-JP  =  0  ,  (2-1) 

where  k0  is  a  reference  wavenumber  and  n(r,z)  is  the  index  of  refraction.  We 
substitute 


,  "/yfy''?-'///"/. 


MEDIUM  1  //////7/7//7/ 


"ui  ,lJl  ,cl 

y//////////y///7/. 


Y///////7 
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P(r,z)  =  u(r,z)H^  (kQ,r)  ,  (2-2) 

into  equation  (2-1).  The  result  obtained  is  the  two-way  (outgoing  and  incom¬ 
ing)  farfield  equation  for  the  acoustic  field,  u, 

urr  +  2ik0ur  +  ko(n2  '  +  u2z  =  0  '  (2-3) 

We  seek  boundary  conditions  for  a  horizontal  interface  that  apply  to  equation 
(2-3). 

The  continuity  of  pressure  gives  rise  to  the  first  interface  condition. 

Pi (r , zB)  =  p2(r,zB)  .  (2-4) 

The  continuity  of  the  normal  component  of  particle  velocity  gives  rise  to  the 
second  interface  condition, 


o 


2 


3Pl 


(2-5) 


By  substituting  equation  (2-2)  into  equation  (2-4) ,  we  obtain  the  first  inter¬ 
face  condition  for  the  parabolic  wave  equation. 


ui (r ’ ZB)  =  u2  (r ’ ZB^ 


(2-6) 


Taking  partial  derivatives  of  both  sides  of  equation  (2-2)  with  respect  to  z 
and  substituting  the  results  into  equation  (2-5),  we  obtain  the  second  inter¬ 
face  condition  for  the  parabolic  wave  equation, 


3ul 
J2  3- 


3u2 
’l  3z 


(2-7) 


2.2  TWO-WAY  FARFIELD  EQUATION  WITH  INTERFACE  CONDITIONS 

The  treatment  to  be  discussed  in  this  section  deals  with  a  horizontal 
fluid  boundary.  A  uniform  partition  in  the  z  direction  is  assumed,  with 
iz  =  h.  The  range  increment  Ar  is  denoted  by  k,  and  m  is  an  integer  index 
that  indicates  the  interface  boundary.  Superscripts  indicate  the  range  level 
and  subscripts  indicate  the  depth  level.  We  use  the  convention  that  if  both 
the  superscript  and  subscript  are  dropped,  the  field  denotes  the  field  at  (nk, 
mh) ,  l  .e.  ,  u  =  u£. 

Consider  the  case  shown  in  figure  1,  where  ulj  denotes  the  field  at 

iti  n+ 1 

z  =  m(Az)  =  mh,  r  =  n(Ar)  =  nk .  Our  main  effort  is  to  determine  u^  such 
that  the  conditions  imposed  by  equations  (2-b)  and  (.2-7)  are  satisfied. 


j> 
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Medium  1 

In  medium  1,  the  field  must  satisfy  equation  (2-3),  i.e., 

(u1)rr  +  2ikQ(u1)r  +  k^(n2  -  l)Uj  +  (Uj)22  =  0  (2-8) 

Using  the  first  three  terms  of  a  Taylor  series  expansion  for  ujj^  upon  u^  and 
solving  for  u,2,  we  find 


(ui  -  ^-i)  +  H -x  TT 


Substituting  equation  (2-9)  into  equation  (2-8)  and  simplifying,  we  obtain 


‘1  h  "  “1 


2  hr2 


+  2ik, 


^  ♦  k2(n2  -  l)Ul  ♦  [-  ^•(u1  -  u^.j]  .  (2- 


Medium  2 

Similarly,  in  medium  2,  the  field  must  satisfy  equation  (2-3),  i.e., 


(u2^rr  +  2ik0(u2)r  +  k0(n2  '  1)u2  +  (u2}zz  =  0 


(2-11) 


Using  the  first  three  terms  of  a  Taylor  series  expansion  for  ujjj+1  upon  u^  and 
solving  for  u.,-,,  we  find 


a-,2  '  '  h 


2  3u2 


h  3  z 


(2-12) 


Substituting  equation  (2-12)  into  equation  (2-11)  and  simplifying,  we  obtain 

lUr,  I  f3  U,  SU^  p 

TT  =  T  7T  +  2ikn  TT  +  kn(n?  "  J)  +  -  u?)  \\  ■  (2‘15) 


In  view  of  the  first  interface  condition,  equation  (2-6),  we  require  that 
u,  =  u.  =  u  in  equations  (2-10)  and  (2-13)  .  Then,  we  multiply  both  sides  of 
equation  (2-10)  by  o2  and  both  sides  of  equation  (2-13)  by  p  and  equate  the 
results,  using  the  second  interface  condition,  equation  (2-7).  After  simpli¬ 
fying,  we  obtain  the  following  parabolic  equation: 


urr  +  2ikour  +  ~ 


♦  — - - (k^  +  —  kf\u  -  k2U 

P1  +  °2\2  P1  7 

2 _ /  P1  *  P2  ^  \ 

p 2)  \  m” 1  p2  u  +  p2Um+1) 


(2-14) 


=  0  . 


.  '  -  *  .  ’  «  *  .  •  -  •  .  •  •  r  .  *  -  •  .  .  .  *.  -  .  *>  ■  v  'i. 
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Equation  (2-14)  is  the  parabolic  wave  equation  satisfied  by  the  field  u 
on  the  interface  Zg.  Note  that  CjCr.z)  =  c0(r,z)  implies  that  kj  =  k2.  If, 
in  addition,  p1  =  p2,  equation  (2-14)  reduces  to 


rr 


2ik  u 
0  r 


k^(n2  -  1 ) u  + 


(um+ i  -  US)  -  (US  -  um- i ) 


=  0 


(2-15) 


Regarding 


h-2 


un) 

(un  -  un  ,)1 

m  / 

l  m  m-i/J 

as  the  central  finite-difference  operator  for  u,,,  equation  (2-15)  becomes 
exactly  the  parabolic  wave  equation,  equation  (2-3),  in  one  medium.  Next, 
define 


and 


and  rewrite  equation  (2-14)  as 

u  ♦  2ik„u  +  Gu  =  0  . 
rr  Or 


(2-16) 


(2-17) 


(2-18) 


Equation  (2-18)  is  the  two-way  farfield  equation  with  interface  conditions 
for  the  acoustic  field. 


2.5  IMPROVED  PARABOLIC  EQUATION 

The  associated  one-way  outgoing  wave  equation  is  obtained  by  using  the 
quadratic  equation  formula  to  formally  solve  equation  (2-18)  for  lu/3r, 

3u/3r  =  +  i  { VM  +  6  -  kju  .  (2-19) 

U  J 

In  equation  (2-19)  ,  the  plus  sign  has  been  chosen  so  that  outgoing  waves  are 
represented.  For  a  range- independent  medium,  equation  (2-19)  is  exact,  that 
is,  a  solution  of  equation  (2-19)  is  also  a  solution  of  equation  (2-3).  Now, 
write 

Q  =  k2  ♦  G  .  (2-20) 

Equation  (2-19)  can  now  be  written  as 

W'<r  =  i(VQ  -  k0)u  .  (2-21) 
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Now  write 


■) 


define 


and  rewrite 


« ■  k5  cf  * 1 


q  =  A.  i, 
4  1-  2 
k0 


Q  =  (q  +  1)  . 


(2-22) 


(2-2 


(2-2 


Using  a  rational  function  approximation  for  the  square  root  operator 


and  substituting  equation  (2-25)  into  equation  (2-21) ,  we  find  that 


3u  . , 

=  ik. 
3r  0 


Rewriting  equation  (2-26) ,  we  have 


,  3u 


•  *)■  • 


(2-2 


(C  ♦  Dq)^  u  =  i k n  [  (A  -  C)  +  (B  -  D)q]u 


(2-2 


Assuming  Q  is  a  constant  over  the  interval  r11  and  rn+1  and  using  the 
Crank-Nicolson  scheme  to  solve  equation  (2-21),  we  find  that 


,n+ 1 


(2-2 


-  u"  =  i(VQ  ' 

where  Ar  =  rn+1  -  rn.  The  resulting  implicit  equation  is 

[l  -  i(  VQ  -  k0)^]un+l  =  [l  ♦  i(  VQ  -  k0)y]un  .  (2-2 

Using  the  definition  for  Q  along  with  the  rational  approximation  for  Q, 
we  obtain 


/  A  +  B  X- 
1  k. 


1  -  i  k  , 


-  1  ^ 


C  ♦  11 


,n+l  - 


/a  + 


1  -  ik. 


C  +  D  —r 
k  0 


'_!L  .  jl  ir 


(2-2 


Substituting  equation  (2-l"j  into  equation  (2-30)  and  simplifying,  we 
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°1  (2  p2  A 

D  - (n2  +  —  nt) 

Lpi +  p2\2  pi  7 


1  -  ikJA  -  C  +  (B  -  D) 


pl  /. 

- I ni 

[pl  +  p2^ 


ni)  -  1 


T  *  pf1'  ’  iko‘“  -  D)fH="n 


=  c  *  D[ *  ^7 "■) 

f  j  -  ‘])t  *  *  ‘ 


—  nj j  -  1  +  ikQ  A  -  C  ♦  (B  -  D) 


—A - 6 

1  + 


(2-31 


+  ik0(B  -  D)^]tz,  un 


Now,  define 


S1  /,  a2  A  , 

n  ■  77V—  p  1 

and  substitute  equation  (2-15)  into  equation  (2-31).  We  now  have 
um+I(C  *  Un  -  ik0t(A  -  C)  -  (B  -  D) n )~~}  +  ^[D  -  ik0(B  -  D)~j 

2?  !  /  °i  +  Pi  \ 

•  - 2 - [un+ 1  -  — - -  un+1  +  —  un+1) 

(dt)2(Dl  +  D2)  \m'1  oz  m  P2  m+1/ 


(2-32 


un(c  +  Dn  +  ikQ  [  (A  -  C)  ♦  (B  -  D)nJ^j  +  ^[d  +  ikQ(B  -  D)^] 


(2-53 


u=r(p 


Next,  define 


2e,  /  +  p2  pi  \ 

— - lu11  -  — - -  Un  +  —  un  J  . 

+  p2) \m-1  p2  m  p2  m+y 


ik„  dr 
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Only  simple  modifications  were  required  to  incorporate  the  improved  PE 
and  horizontal  interface  conditions  into  the  IFD  model. ® 


2.4  WIDE-ANGLE  CAPABILITY 

The  capability  to  handle  wider  propagation  angles  results  from  judicious 
selection  of  the  coefficients,  A,  B,  C,  and  D  in  the  rational  function  approx¬ 
imation  for  the  square  root  operator  in  equation  (2-25). 

If  we  select  A  =  1,  B  =  1/2,  C  =  1,  and  D  =  0,  equation  (2-41)  reduces  to 
a  solution  of  the  standard  Tappert  parabolic  equation.  Tappert's  equation  is 
limited  to  propagation  angles  of  less  than  approximately  15  deg. 

Claerbout ' s*4 » 15  one-way  equation  uses  A  =  1,  B  =  3/4,  C  =  1,  and  D  =  1/4. 
Claerbout's  coefficients  extend  the  propagation  angle  to  approximately  40  deg. 
Claerbout 's  approximation  is  equivalent  to  a  quadratic  expansion  of  VQ, 
whereas  Tappert's  approximation  is  a  linear  expansion. 


3.  COMPUTER  IMPLEMENTATION 


The  model  that  implements  the  implicit  finite-difference  (IFD)  formula, 
equation  (2-41),  has  been  written  in  FORTRAN  using  single-precision  complex 
arithmetic  and  has  been  installed  on  a  VAX-11/780  digital  computer.  For  a 
detailed  description  of  the  IFD  model,  the  reader  is  referred  to  reference  9. 

The  next  two  subsections  describe  input  and  output  formats  in  detail.  The 
only  difference  in  formats  between  this  and  the  original  IFD  model  is  in  the 
insertion  of  card  2A  into  the  input  runstream.  Test  examples  showing  a  variety 
of  sample  input  runstreams  and  user-written  subroutines  can  be  found  in  refer¬ 
ence  9.  To  rerun  those  test  examples,  card  2A  must  be  inserted  in  the  run- 
stream  . 


3.1  INPUT  FORMAT 

Prior  to  executing  the  IFD  model,  input  card  images  containing  problem 
parameters  must  be  stored  in  file  IFD. IN.  File  IFD. IN  is  assigned  to  FORTRAN 
unit  number  NIU  in  the  main  program.  If  the  user  prefers  to  input  problem 
parameters  on  cards,  then  parameter  NIU  should  be  equated  to  the  card  reader 
unit  number,  and  the  statement  that  assigns  file  IFD. IN  should  be  removed  from 
the  main  program.  In  either  case,  the  input  runstream  is  prepared  in  free  for¬ 
mat  as  follows. 


9 


TR  6905 


CARD  CONTENTS 


1  FRQ,  ZS,  CO,  ISF,  RA,  ZA,  N,  IHNK,  ITYPEB,  ITYPES 

where 

FRQ  =  frequency  (Hz) . 

ZS  =  source  depth  (m) . 

CO  =  reference  sound  speed  (m/s) . 

If  CO  =  0.0,  CO  is  set  to  the  average  sound  speed  in  the  first 
layer . 

ISF  =  starting  field  flag. 

0  =  Gaussian  starting  field  is  generated. 

1  =  user  prepares  starting  field.  See  subroutine  UFIELD. 

If  ISF  =  0,  RA  is  set  to  zero. 

RA  =  horizontal  range  from  source  to  starting  field  (m) . 

If  ISF  =  0,  RA  is  set  to  0.0. 

ZA  =  depth  of  starting  field  at  range  RA  (m) .  If  ZA  =  0.0,  ZA  is 
set  to  the  maximum  depth  of  the  bottommost  layer  in  the  first 
profile. 

If  ITYPEB  =  2  or  3  and  ZA  =  0.0,  ZA  is  set  to  (4/3)*  maximum 
depth  of  the  bottommost  layer.  If  ITYPEB  =  2  or  3  and  ZA  /  0.0 
the  artificial  bottom  layer  is  extended  to  ZA  meters  provided 
that  ZA  is  greater  than  or  equal  to  the  maximum  depth  of  the 
bottommost  layer  in  the  first  profile. 

N  =  number  of  equispaced  receivers  in  the  starting  field.  If 

N  =  0,  N  is  set  so  that  the  receiver  depth  increment  is  less 
than  or  equal  to  1/4  wavelength.  If  N  is  greater  than  MXN,  N 
is  set  to  MXN.  See  parameter  MXN. 

IHNK  =  Hankel  function  flag.  If  IHNK  =  0,  don't  use  Hankel  function. 

If  IHNK  =  1,  divide  the  starting  field  by  the  Hankel  function, 
then  multiply  the  solution  field  by  the  Hankel  function  before 
computing  propagation  loss.  If  starting  field  is  Gaussian, 

IHNK  should  be  set  to  0.  If  starting  field  is  elliptic,  IHNK 
should  be  set  to  1 . 


ITYPEB  =  type  of  bottom. 

=  0,  homogeneous  Neumann  boundary  condition;  program  supplies 
bottom  condition. 


CONTENTS 


=  1,  user  supplies  bottom  condition.  See  subroutine  BCON. 

=  2,  artificial  absorbing  layer  introduced;  bottom  of  layer 
follows  contour  of  water-bottom  interface. 

=  3,  artificial  absorbing  layer  introduced;  bottom  of  layer 
kept  flat. 

ITYPES  =  type  of  surface. 

=  0,  pressure  release;  SCON  sets  SURY  and  SURX  =  0.0. 

/  0,  user  inserts  code  in  SCON  to  compute  SURY  and  SURX. 

RMAX,  DR,  WDR,  WDZ,  PDR,  PDZ,  ISFLD,  ISVP,  IBOT 

RMAX  =  maximum  range  of  solution  (m) . 

DR  =  range  step  for  marching  solution  Cm) . 

If  DR  =  0,  DR  is  set  to  1/2  wavelength. 

WDR  =  range  step  (rounded  to  nearest  DR)  at  which  solution  is  writ¬ 
ten  on  disk  (m) . 

WDZ  =  depth  increment  (rounded  to  nearest  DZ)  at  which  solution  is 
written  on  disk  (m) . 

PDR  =  range  step  (rounded  to  nearest  DR)  at  which  solution  is 
printed  (m) . 

PDZ  =  depth  increment  (rounded  to  nearest  DZ)  at  which  solution  is 
printed  (m) . 

ISFLD  =  0,  don't  print  starting  field. 

=  1,  print  starting  field. 

ISVP  =  0,  don't  print  sound-speed  profile. 

=  1,  print  sound-speed  profile. 

IBOT  =  0,  don't  print  bottom  depths. 

=  1,  print  bottom  depths. 

A,  B,  C,  D 

TAPPERT:  A  =  1 ,  B  =  0 . 5 ,  C  =  1 ,  and  D  =  0 . 
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CONTENTS 


CLAERBOUT:  A  =  1.0,  B  =  0.75,  C  =  1.0,  and  D  =  0.25. 


Rl,  Z1 
R2 ,  12 


R3 ,  Z3 


Bottom  profile.  Range  and  depth  of  water  (m) ;  maximum 
number  of  depths  =  100  (see  parameter  MXTRK) . 


N+l  -1,  -1 


N+2  RSVP 


N+3  KSVP 


N+4  NLYR 


Required.  Marks  end  of  bottom  profile. 


N+5  ZLYR(I) ,  RHO(I),  BETA(I) 

N+6  ZSVP(l),  CSVP(l) 

N+7  ZSVP (2) ,  CSVP(2) 


Repeat  for  each 
layer.  1=1,  NLYR. 


Repeat  for  each 
profile. 


N+M  ZSVP(J),  CSVP(J) 


where 


=  range  of  SVP  (m) . 

=  SVP  flag. 

=  0,  profile  is  in  runstream. 

/  0,  profile  (cards  N+4  through  N+M)  is  supplied  by  user-writ¬ 
ten  subroutine  USVP.  KSVP  may  be  used  in  computed  GOTO 
statement  to  transfer  control  in  user  subroutine  USVP. 


NLYR  =  number  of  layers,  if  ITYPEB  =2  or  3,  the  program  inserts  an 
artificial  absorbing  layer  and  then  increments  NLYR  by  1. 
Maximum  NLYR  =  100  (see  parameter  MXLYR) . 

ZLYR(I)  =  maximum  depth  of  layer  I  in  profile  (m) . 
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CARD 


CONTENTS 


RHO(I) 

BETA(I) 

ZSVP 

ZSVP(l) 

ZSVP(J) 

CSVP 

CSVP(l) 

CSVP(J) 


density  in  layer  I  (g/cm3). 
attenuation  in  layer  I  (dB/wavelength) . 

SVP  depth  (m) .  Maximum  number  of  sound-speed  profile  values 
=  100  (see  parameter  MXSVP) . 

depth  to  top  of  layer  I  (m) . 

depth  to  bottom  of  layer  I  (m) . 

SVP  speed  (m/s).  Maximum  number  of  sound-speed  profile  val¬ 
ues  =  100  (see  parameter  MXSVP). 

speed  of  sound  at  top  of  layer  I  (m/s) . 

speed  of  sound  at  bottom  of  layer  I  (m/s). 


3 . 2  OUTPUT  FORMAT 

Output  from  the  IFD  model  is  written  on  disk  in  file  IFD.OUT.  File 
IFD.OUT  is  assigned  to  FORTRAN  unit  number  NOU  in  the  main  program.  The  data 
written  in  IFD.OUT  are  unformatted  and  are  written  with  FORTRAN  WRITE  state¬ 
ments  as  follows. 

WRITE (NOU) FRQ , ZS ,C0 , ISF , RA, ZA,N , IHNK, ITYPEB, ITYPES ,RMAX,DR,WDR,DZ, 

NLYR, ZLYR,RHO,BETA. 

WRITE(NOU)NZ,RA,WDZ, (U (I ) , I=IWZ,N, IWZ) . 

The  first  WRITE  statement  is  executed  only  once  and  writes  the  value  of 
each  of  the  following  parameters  at  the  start  of  the  problem. 

FRQ  =  frequency  (Hz) . 

ZS  =  source  depth  (m) . 

CO  =  reference  sound  speed  (m/s) . 

ISF  =  starting  field  flag. 

=  0,  Gaussian. 

=  1,  user. 

RA  =  horizontal  range  from  source  to  starting  field  (m) . 

ZA  =  depth  of  starting  field  (m) . 

N  =  number  of  equispaced  receivers  in  the  starting  field. 


l.i  ^  ^  ^  ^  . 


♦  *  v *  * 

*  -V 
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IHNK  =  Hankel  function  flag. 

=  0,  Hankel  function  not  used. 

=  1,  starting  field  was  divided  by  Hankel  function. 

ITYPEB  =  type  of  bottom. 

=  0,  rigid. 

=  1,  user  bottom. 

=  2,  artificial  absorbing  layer,  follows  bottom  contour. 

=  3,  artificial  absorbing  layer,  bottom  kept  flat. 

ITYPES  =  type  of  surface. 

=  0,  pressure  release. 

+  0,  user  supplied. 

RMAX  =  maximum  range  of  solution  (m) . 

DR  =  range  step  for  marching  solution  (m) . 

WDR  =  range  step  at  which  solution  is  written  on  disk  (m) . 

DZ  =  depth  increment  of  receivers  (m) . 

NLYR  =  number  of  layers. 

ZLYR  =  array  containing  depth  of  each  layer  (m) . 

RHO  =  array  containing  density  in  each  layer  (g/cm3). 

BETA  =  array  containing  attenuation  in  each  layer  (dB/wavelength) . 

The  second  WRITE  statement  is  executed  at  each  write-range  increment,  WDR 
The  data  written  are  as  follows: 

NZ  =  number  of  equispaced  receivers  in  the  solution  field. 

RA  =  horizontal  range  from  source  to  solution  (m) . 

WDZ  =  depth  increment  at  which  solution  is  written  on  disk  (m) . 

U  =  array  that  contains  the  complex  field  at  range  RA. 

If  IHNK  =  1,  then  the  contents  of  U  must  be  multiplied  by  the 
Hankel  function  before  computing  propagation  loss. 

IWZ  =  index  increment  of  receiver  solutions  to  be  written  on  disk. 


The  following  READ  statement  may  be  used  to  read  the  solution  field: 
READ (unit)  NZ, RA,WDZ, (U(I) 1=1 ,NZ) . 


4.  TEST  PROBLEM 


The  new  IFD  model  with  wide-angle  capability  is  tested  with  test  case  3B, 
extracted  from  reference  16.  The  environment  for  this  case  is  range-independ¬ 
ent  and  consists  of  an  isovelocity  water  column  over  an  isovelocity  half-space 
bottom.  Parameters  for  the  problem  are  given  below: 

Frequency  =  250  Hz, 

Water  depth  =  100  m, 

Sound  speed  in  water  =  1500  m/s. 

Density  in  water  =  1.0  g/cm3. 

Density  in  bottom  =  1.2  g/cm3. 

Attenuation  in  water  =  0, 

Bottom  attenuation  =  0.5  dB/X, 

Bottom  sound  speed  =  1590  m/s. 

Maximum  range  =  10  km. 

Source  depth  =  99.5  m. 

Receiver  depth  =  99.5  m,  and 
Modes  =  11. 

With  the  source  and  receiver  just  off  the  bottom,  the  higher  modes  are 
more  strongly  excited,  and  all  11  modes  are  required  for  an  accurate  solution 
in  the  5  to  10  km  range  interval.  This,  in  turn,  requires  a  half-beamwidth 
capability  of  18.5  deg. 

Comparisons  with  normal-mode  results  indicate  that  the  fast-field  program 
(FFP)  solution  shown  in  figure  2  is  correct. H>17 

Two  IFD  solutions  to  test  case  3B  were  obtained,  one  using  Claerbout's 
coefficients  A  =  1,  B  =  3/4,  C  =  1,  and  D  =  1/4  (wide-angle  PE),  and  the  other 
using  the  coefficients  A  =  1,  B  =  1/2,  C  =  1,  and  D  =  0  (standard  Tappert  PE). 
In  each  case,  a  Gaussian  pressure  distribution  and  an  artificial  boundary 
extended  to  250  m  in  depth  were  assumed.  As  shown  in  figure  3,  the  IFD  solu¬ 
tion  obtained  using  Claerbout's  coefficients  differs  slightly  in  level  but,  in 
general,  is  in  excellent  agreement  with  the  FFP  solution.  The  average  level 
of  the  IFD  solution  to  the  standard  PE  is  in  good  agreement  with  the  FFP  but 
the  pattern  of  the  solution  is  poor. 


The  input  runstream  for  test  case  3B  is  listed  below: 


250  99.5  0  0  0  250  1000  030 

10000  1  2  99.5  1000  20  0  0  0 

1  .75  1  .25 
0  100 

10001  100 

-1,  -1 

0 

0 

2 

100  1  0 
0  1500 
100  1500 
200  1.2  0.5 
100  1590 
200  1590 


5.  CONCLUSIONS 

The  mathematical  formulation  of  a  recent  improvement  that  extends  the 
modeling  of  propagation  angles  from  less  than  approximately  15  deg  up  to 
approximately  40  deg  and  the  modifications  that  incorporate  this  wide-angle 
capability  into  the  1FD  model  were  presented. 

IFD  calculations  with  and  without  the  wide-angle  capability  were  com¬ 
pared  with  normal-mode  and  FFP  solutions  for  an  isovelocity  range- independent 
test  case.  IFD  calculations  obtained  with  the  wide-angle  improvement  are  in 
excellent  agreement  with  the  normal-mode  and  FFP  solutions,  whereas  IFD  calcu¬ 
lations  obtained  without  the  improvement  are  subject  to  phase  errors. 

As  demonstrated  in  this  report,  the  capability  to  handle  wide-angle 
propagation  was  incorporated  into  the  IFD  model  quite  easily.  Future  capa¬ 
bilities  to  be  built  into  the  model  include  multiple  irregular  interfaces, 
automatic  step-size  determination,  high-frequency  propagation,  shear  waves, 
and  others. 
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